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Abstract 

We consider a nonlinear field equation which can be derived from a binomial lattice as 
a continuous limit. This equation, containing a perturbative friction-like term and a free 
parameter 7, reproduces the Toda case (in absence of the friction- like term) and other equa- 
tions of physical interest, by choosing particular values of 7. We apply the symmetry and 
the approximate symmetry approach, and the prolongation technique. Our main purpose 
is to check the limits of validity of different analytical methods in the study of nonlinear 
field equations. We show that the equation under investigation with the friction-like term is 
characterized by a finite-dimensional Lie algebra admitting a realization in terms of boson 
annhilation and creation operators. In absence of the friction-like term, the equation is 
linearized and connected with equations of the Bessel type. Examples of exact solutions 
are displayed, and the algebraic structure of the equation is discussed. 
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1 Introduction 



In this paper we study the class of nonlinear wave equations 
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where u = u(x,t), subscripts denote partial derivatives, e is a free parameter, and 
7 is a constant which takes all (real) values except zero and one. 

For e = 0, Eq. (|l.l|) can be considered as the continuous limit of a uniform one- 
dimensional nonlinear lattice of N particles interacting through the nearest-neighbour 
potential ffl] 
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where a n and b n are constants of the n-th nonlinear spring. 
The equation of motion of the chain reads 

M n y n — — [(j>'(r n+1 ) — (f>'{r n )) , 



;i.3) 



where <f>' = M n is the mass of the n-particle and r n = y n — y n+ \, y n being 
the displacement of the n-particle from its the equilibrium position. We shall call 
"binomial lattice" the chain governed by the potential ( |1.2|) . If the binomial lattice 
is uniform, i.e. a n = a, b n = b and M n = M for any n, the equation of motion ( |1.3| ) 
takes the form 
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For 7 — > oo, 7 = 2, 7 = 3 and 7 = —1, Eq. (|1.4|) reproduces the Toda @, the 
harmonic, the Fermi-Pasta-Ulam and a Coulomb-like uniform lattice, respectively. 
Since Eq. ( |1.4|) models interesting physical situations, its analysis may be important 
mostly from the point of view of a unifying theory of lattice systems described by the 
potential (|L2|) . 

This difficult task is made simpler by starting from the investigation of the con- 
tinuous limit of Eq. (|1.4| ). Then, Eq. (|1.4j) becomes 
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after the rescaling br — > u and Jjzt-tt. 
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For 7 — » oo, Eq. (|1.5|) coincides with the (l+l)-dimensional continuous Toda equa- 
tion g 

««=(e tt )„, (1.6) 
which is known also as heavenly equation and plays a crucial role in General Relativity 
i, 1- 

Equation ( |1 . 1|) represents a modified version of Eq. (|1.5| ). It contains a friction-like 
term which can be considered as a small perturbation. This extension of Eq.(1.5) 
allows us to check the validity of different analytical methods usually employed to 
handle integrable equations. 

Precisely, here we apply to Eq. (|l.l| ) the symmetry |7[ and the approximate sym- 
metry approach ||, and the prolongation procedure 0. 

The main results achieved in this paper are the following. In Sec. 2 we find the 
approximate symmetry generators of Eq. ( |1 . 1| ) in the case in which e is a perturbative 
parameter. The commutation relations among these operators do not define an exact 
finite-dimensional algebra. (Exact algebras arise only for e = and when e 7^ is 
not considered as a small quantity). Notwithstanding, an exact finite-dimensional Lie 
algebra can be constructed by introducing certain auxiliary operators. A realization of 
this algebra, which characterizes the approximate symmetries of Eq. ( |1.1| ), is obtained 
in terms of boson annihilation and creation operators. 

In Sec. 3, the symmetry approach is exploited to determine special solutions of 
Eq. ( |1 . 1|) and Eq.(1.5). In this context, the study of the symmetry corresponding to 
the generator X°4 (see (2.5c)) is of particular interest. In fact, X°4 is a scale symmetry 
in the variables x, u together with a translation of along u. For 7 — > 00, i.e. when 
Eq.(l.l) coincides with the continuous Toda equation, the scale symmetry is broken. 
The reduced equation related to X° 4 resembles the one-dimensional Liouville equation 
(see (3.4)) containing a square derivative which vanishes in the limit 7 — > 00. In the 
case e = 0, Eq.(3.4) is exactly solved for some values of the parameter 7. For example, 
if 7 = 3 (a choice associated with the Fermi- Past a- Ulam potential), we get a solution 
to Eq.(1.5) via Eq.(3.4) expressed in terms of the Weierstrass function. 

For e ^ 0, we provide an implicit solution to Eq.(l.l) by using the symmetry 
variable V = vd x + d t , where v is a constant. This solution (see (3.42)) involves a 
generalization of the exponential-integral function. It is notable that for 7 — > 00, this 
solution reproduces just that corresponding to the continuous Toda case [[K] 



Sec. 4 is devoted to a study of Eq.(l.l) within the prolongation scheme 0. This 
method is particularly convenient because it allows us to analyze Eq.(l.l) for e = 
and e 7^ by an algebraic point of view. We show that the (differential) prolongation 
equations for Eq.(1.5) afford a class of solutions connected with the Lie algebra of 
the Euclidean group E2 in the plane. This algebra leads to a linearized version of 
Eq.(1.5), in the sense that a mapping is established between Eq.(1.5) and the linear 
wave equation (4.6) written in a pseudopotential variable. For 7 — > 00, Eq.(4.6) 
coincides with a two-dimensional form of a linear wave equation that occurs in a 
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generalized Gibbson-Hawking ansatz On the other hand, the prolongation 



equations for Eq.(l.l) can be solved in terms of power series expansions whose coef- 
ficients (in the variable z = (l + ) depend on the pseudopotential components 
and, for a certain infinite set of values of 7, obey a finite-dimensional Lie algebra. A 
representation of this algebra indicates the existence of a possible link between the 
prolongation method and the symmetry approach. Finally, in Sec. 5, some comments 
are reperted. 



2 The approximate symmetry approach 

In the case in which e is a perturbative parameter, Eq.(l.l) can be handled profitably 
by means of the method devised by Baikov, Gazizov and Ibragimov some years ago 
||. This method enables one to construct approximate symmetries which are stable 
for small perturbations of the differential equation under investigation. 

In order to obtain the approximate symmetries of Eq. (1.1) (e 7^ 0), we have to 
introduce the quantities 



^0 = u tt - 
Fi = u t , 



[1 + u/j) 
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into the equations (3.16) and (3.17) of ||. 

For brevity, we shall omit any calculation and report below the results only. 
The approximate symmetry generator turns out to be 



X 



h- 2 * , 

cit + c 2 + e — ~t;Ci— + kit + k 2 
V7 + 2 2 



+ 



7-2 V 7/ 7 



7 



d t + [(ci + c 2 ) x + c 4 + e ((hi + fc 3 ) x + fc 4 )] d x 
7-2 



7 + 2 



1 + H 

7, 



(2.2) 



with ci, C2,... and k\, k 2 ,... arbitrary constants. 
From (TT) we get 



X x = X 1 ° + e 
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where 



7 + 2 2 7 + 2 
Xl = td t + xd X} 
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(2.5e) 
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The operators X° are the (exact) symmetry generators of Eq. (1.1) 

for e = 0, while X\, X 2 , X 3 , X 4 are the approximate symmetry generators of Eq. 
(1.1) for e 7^ 0. The exact symmetry generators of Eq.(l.l) for e ^ are X®, X° and 
X 3 . The operators X®, X$, X®, X® are inessential. 

The approximate symmetry generators X\, X2, X 3 , X 4 satisfy the commutation 
relations 



[*2,*i] 

[x 2 ,x 4 ] 

[Xi,X 3 ] 



X 2 + e 
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td t 
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[X 2 ,X 3 ]=0, [X h X 4 ] 
X3, [X 3 , X 4 ] = x 3 



(2.6) 

(2.7) 
(2.8) 



The commutation rules ( |2.6|) - (|2.8|) define an exact finite-dimensional Lie algebra 
only for e = 0. However, they can be exploited to build up a realization of an exact 
Lie algebra, which holds for e 7^ 0, with the help of the "auxiliary" operators 
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In doing so, we find 
[X 2 ,X l ] 

[x 2 ,x 4 ] 

[XuX 3 ] 
[Xi,Z] 

[*2,Y] 

[X 3 , Y] 



7 + 2 2 



7 + 2 \ 7, 



d u . 



X 2 + eZ, 

[X 2 ,X 3 ] = [X 1; X 4 ] =0, 

— ^3> [X 3 ,X 4 ] = X 3 , 

-e^Y, [X U Y)=Y, 

Z, [X 2 ,Z] = 1^X 2 , [Z,Y] 
[X 3 ,Z] = [X 4 ,Y] = [X 4 ,Z]=0. 



7 



-Y. 



(2.9a) 
(2.9b) 

(2.10a) 
(2.10b) 
(2.10c) 

(2.10d) 

(2.10e) 
(2.10f) 
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Hereafter, the symbols Y, Z, Xj (j = 1, ..,4) will mean both the abstract ele- 
ments and the corresponding realizations (see (|2.5|) , ( |2.9|) and (|2.10|) ) of the finite- 
dimensional Lie algebra ( |2. 10| ) ■ 

At this stage, we remark that the subalgebra ( |2.10e ) is isomorphic to the s/(2, R) 
algebra 

\Z', T] = 2S, [T, S] = 2Z', [S, Z'\ = -2T, 



(2.11) 



where 



T^2l±| (Y + X 2 ), 

Z> = 2^ Z. 

7-2 

It can be proved straightforwardly that 
Proposition 1 The Casimir operator 

C = T 2_ S 2_ z ,2 
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(2.12) 
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of the Lie algebra ( \2.11\ ), commutes with all the elements Y, Z, Xj ( j = 1, ..,4) of the 
Lie algebra ( \2.1Q ). 

Furthermore, the following Proposition holds: 



Proposition 2 The Lie algebra ( \2. 1 Q ) admits a realization in terms of boson anni- 
hilation and creation operators. 

This can be seen by setting 

a+ = t, a\ = u, aj = x, a x = d t , a 2 = d u , a 3 = d x , 



to give 
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(j) k = 1,2, 3), (2.15) 
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(2.16) 
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and 
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where = — 8a 2 (2a 2 + 1) is the Casimir invariant relative to the continuous 
Toda case 7 — > 00, and the operator C (7), defined by 



C( 7 ) 
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tends to zero as 7 — ► 00. 

Although the full role of the closed algebra (|2.10| ) has to be better understood, it 
is noteworthy that the auxiliary operators Z and Y (see ( |2.9| )) can be interpreted as 
symmetry variables of the equations 



7 



1 + 



u 



u., 



7, 



u 



7, 



1 + 



II 



1-7 



7, 



3 / u 
u xx + eu x - — [1 + - 
27 V 7, 



-1 



ui = -2 1 + 



7, 



1-7' 



1 + 



u 



7, 



(2.19) 



(2.20) 



respectively. 

Equations ( |2.19| ) and Q2.20| ) arise formally from Eq.(l.l) via the transformations 
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(2.22) 



From Z we obtain the invariants x' = x, r]'(x') = r](x), with 
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(2.23) 



which gives the reduced equation 
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(2.24) 



from substitution into (|2.19|) . We notice that Eq.( [2.24j ) coincides, formally, with 
Eq.( p.3|) for x — > t (see Sec. 3). Therefore, putting 77 = e^ 7 , Eq. ( |2.24| ) becomes 
Eq.( [3.4|) where t is replaced by x. 
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On the other hand, the operator Y yields the invariants x' 
with 



r(x) 



x . 



r'(x') 




r[x) 



(2.25) 



which leads to the ordinary differential equation of the modified Liouville type (see 

ID) 



+ ev x + 
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oo, in the reduced 



7-2 x 7-2 

where v = lnr. 

We point out that in the continuous Toda case, i.e. for 7 
equations of the one-dimensional modified Liouville type [ (p.4j) , (|3.4j) written in the 
variable x and ( [2.26| )1, the terms Wj 2 , W%, and v\ disappear. This is due to the fact 
that for 7 — > 00, the coefficient in front of ud u in X° is vanishing, so that X% is 
not longer a scale symmetry in the variables x, u (together with a translation of 
-^2 along u). Thus, the presence of square first derivatives in the reduced equations 
generated by X% reflects the existence of a scale symmetry, which is broken in the 
limit 7 — > 00. 



3 Explicit solutions 

Here we shall display some significant examples of exact solutions to Eq. (|1.1|) by using 
the symmetry approach. To this aim, first let us deal with the generator X° which 
appears in both the cases e = and e / 0. The group transformations involved by 
Xl are 



x 

t' = t 



e x x, 



u 
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1 + -)e^- 2 
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where A is the group parameter. A set of basis invariants corresponding to 
given by 



1 + — ) 2 x 

7 



1 + -) 2 X) 
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h = 
h = t' = t. 

Now, by putting p(t) = J 1; from Eq. (|1.1| ) we obtain the reduced equation 

2 2-7 2-7 7-I 
Pt + PPtt + epp* = • 

7 7 7 



(3.1a) 
(3.1b) 

(3.1c) 
0) is 

(3.2a) 
(3.2b) 

(3.3) 
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Through the change of variable p = e , Eq.( |3.3| ) becomes 



W tt + e W t + 
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(3.4) 



which is a kind of modified one-dimensional Liouville equation. For e = 0, Eq.(|3] 
can be solved exactly. This occurs via the position 



W 
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hie. 



which transforms Eq.( |3.4| ) into 
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(2- 7 ) 2 



(3.5) 
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From (|3.6[) one easily finds 



(2 - 7) 2 7 

where c is a constant of integration. 
Putting 

y 2 = i + a en 



with 



Eq.Q3.7D provides 



4(7-1) 
07(7 - 2) 



2 - 



l— 7 i 

dy (y 2 — 1)~ = — a~ \fc t + const, 



(3.7) 

(3.8) 
(3.9) 

(3.10) 



with c > 0. 

In some situations (see Case iv), it may be convenient to write formula ( J3.1U ) as 



2 7 1 

(coshz)~~ dz = (— 1)~ — a~ y/c t + const, 



(3.11) 



which follows from ( |3.10 ) via the change of variable y = tanh^. 

Once the integral at the left-hand side of ( |3. 10|) has been calculated, one obtains 
the exact solution 

y - 1\" 



u = 7 



X j-2 



(3.12) 



to Eq.flLlD (e = 0) [see (|3T2a|) , (gg), and (|3^) 1. 

At this point, by way of example, we would like to deal with some particular value 
of the parameter 7 which is compatible with an explicit expression of y in terms of t. 
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In what follows we shall deal with the cases i) 7 = |, ii) 7 = 3, and Hi) 7 = —2. 

The last two choices correspond, respectively, to the Fermi-Pasta-Ulam potential ||, 
and to a potential whose nonlinear part, 1/(1 — |) 2 , mimics the inverse square po- 
tential appearing in the treatment of the scattering states in conformally invariant 



Quantum Mechanics p~2 |. 

Case i) For 7=5, from ( |3.10|) we have 



J (y 2 ~ l )dy = ^a 2 Vct, 



(3.13) 



with a — — P, where the constant at the right-hand side of ( |3.10| ) has been put equal 
to zero. 

Eq.( |3.13| ) yields the cubic equation 



y 3 - 3y - -a 2 ^~ct = 0, 
which affords the real solution ([13|], p. 17) 

y = s 1 + s 2 

where 

'3 



si 



s 2 



-a 2 v/ct + \ —a 4 ct 2 - 1 
V 64 



(3.14) 



(3.15) 



(3.16a) 
(3.16b) 



for t > 



3a 2 y/c' 



by 



The remaining solutions of (|3.14| ) are complex conjugate functions, and are given 



yi = ~(s 1 + s 2 ) + i-^-(s 1 -s 2 ), y2 = y\- (3.17) 

Now, by inserting ( 3.17 ) in ( |3.14 ) for 7 = ■= we find the exact solution 



1 



u 



— x-l (sl + s 2 2 + l) - 1 



Case ii) For 7 = 3, Eq.( [3.7|) provides 



3 + C, 



(3.18) 



(3.19) 
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which can be written as 



via the rescaling r 
the equation 



9 2 = 4£ 3 - g 3 , (3.20) 
t, with g 3 = -c/2. Equation flOol) is a special version of 



0? 



4 #3 _ _ ^ 



(3.21) 



which is satisfied by the Weierstrass elliptic function V(r; g2, #3), where g2 and #3 are 
the invariants of V (see J3T3H , p. 918). 

Concerning Eq. ( |3.20| ) we shall distinguish two cases, g 3 = and g 3 7^ 0. For 
g 3 = 0, we have 

0(r) = 7>(r;0,0) = 1 (3.22) 

which coincides with the first term of the series representation of V(t; g2, g^)- 
Then, keeping in mind fl3.8|) and ( |3.14j ), we find 



■it 



For g 3 ^0, Eq.(g20|) gives 



da 



that is to say 

Hence, from (|3.12|) we deduce 

u = 3 (x 2 e 



00 y/Aa 3 — #3 
= 7>(r; 0,(73). 



T, 



^ 2 ^(tjO,^) 



(3.23) 

(3.24) 
(3.25) 

(3.26) 



We notice that for c = —2 (g 3 = 1), the Weierstrass function in Q3.26 ) reduces to 
the equianharmonic case V{r\ 0, 1) (JTB|, p. 652). 
Case iii) If 7 = —2, from ( |3.10| ) we easily obtain 



'cfc 1 



X 



(3.27) 



with a = 

8c 



Case iv) For 7 = —1, the potential (|1.2j ) takes the form <fi ~ (1 — it) -1 — (1 + u), 
which resembles a special case of the generalized Killingbeck potential V = — — + 
Br + Cr 2 |14[|. We remind the reader that a Coulomb potential perturbed by a linear 
term describes the spherical Stark effect in hydrogen [|15|]. Putting in (3.11) 7 = — 1 



and c = — |c|, we arrive at the relation 



• , . 9, .3 
2z + smh 2z = -\c\ 2 t + const. 
4 



(3.28) 
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Curiously enough, by replacing formally the hyperbolic sine by the exponential, 
Eq. fl3.28D becomes an equation of the Schroeder type, which appears in the bootstrap 
model and in renormalization theory, whose analytical structure was investigated by 
Hagedorn and Rafelsky ]TjJ. These authors obtained a solution of the Schroeder 
equation both as a power series expansion and as an integral representation. Thus, 
it should be of interest to try to adopt the same strategy in the study of Eq. (|3.28|) . 

Another nontrivial example of exact solution to Eq.(l.l) can be determined start- 
ing from the linear combination of X3 = X° = d x and X2 = X° = df. 



where v is a (real) constant. The group transformations are x' = x + v A, t' — t + A, 
which provide the invariant £ = x' — vt' — x — vt. Thus, by inserting the variable £ 
into Eq.(l.l), we get the ordinary differential equation 



where Co is a constant of integration. By choosing cq = 0, v 2 = :L ^- and limiting 
ourselves to consider those values of 7 such that 7 = — 17|, Eq. fl3.3lD yields 



V = vd x + d t , 



(3.29) 




(3.30) 



which gives 




(3.31) 




(3.32) 



from which 




) 



(3.33) 



£0 being a constant of integration. 



We notice that for I7I — > 00, Eq. (|3.33|) becomes 




(3.34) 



where Ein(u) denotes the exponential-integral function ( ||17|| , p. 255) 



Ein(u) 



Ei(-u) + lnw + C, 



(3.35) 
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Ei(—u) = — s —r du', and C is the Euler-Mascheroni constant, denned by [[L7J 



C = -ij(l) 



e In t dt, 



(3.36) 



where ip(z) = is the psi (or digamma) function p. 258). 
At this point, let us introduce the function 



"(2+|7l) 



{u') a ~ X du' 



u 



"(2+l7l) 



(u') a ~ du' (3.37) 



T(a,u; | 7 |) = 

= i>,0;| 7 |) 

which holds for < a < 2 + | 7 |. 

We point out that for | 7 | — > oo, r(a,w; |7|) reproduces the incomplete gamma- 
function F(a,x) ( |13| p. 260). 

Now, since the integral at the right-hand site of ( |3.37| ) can be written as [ 17 ] 



[u 



/\a-l 



1+ k 



-(2+I7I) 



Vj 11 

du' = — 2 F 1 {2 + | 7 |, a; a + 1, -— ) 
a \j\ 



u 



a 00 



r(2 + | 7 |) 



^ r(2 + | 7 |+n) ( n , 



n=0 



(a + ri)n\ 



H 



(3.38) 



Eq.( [3.37|) takes the form 



r(a, M ;| 7 |) = r(a,0,| 7 |) 



o 00 



^ r(2 + | 7 |+n) ( i)r / // 



r(2+| 7 |)£i (a + n)n! 



(3.39) 



where T(-) is the gamma-function and 2-^1 is the the Gauss hypergeometric series, 
respectively ( |13|] , p. 556). 



We have (see JT7J, p. 209) 

lim a ^ f (a, it; | 7 |) = 

"af (a,0; | 7 |) - w Q 



lim 



a 00 



a 



^ r(2 + | 7 | + 7i) ( l)n 



r(2 + | 7 |)^ (a + n)n\ 

/ \ -(2+I7I) 



ln M + ln| 7 |-^(2+| 7 |)-C-/ o U i + 



du'. (3.40) 
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On the other hand, from ( |3.37| ) we obtain (see [13| p. 255) 

lim^o f(a,«; |7|) = 
Finally, by comparing (|3.41|) with ( p.40| ), we find 



00 / W 

l 1 + R 



-du'. 



w 



' u'Y (2+h]) 1 hi 

1 + rT -du' = ln^-^(2+| 7 |)-C + 

n\J u u 



ru 1 


1 - 
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Jo u' 







"(2+l7l)' 



that is (see (|3.33|) ) 



r(o,« ; | 7 |) 



1 + — — du' 

v 7 y « 



-(€ - Co) + In ^ - ^(2 + |7l) - C. 
v u 



We remark that lini| 7Hoo [In | 7 | - ^(2 + | 7 |)] = (pf, p. 945). 



Therefore, Eq. ( |3.43| ) becomes 



lim^i^oo f(0,u;| 7 |) = -Ei(-u) = -(£ - f ) - Inn - C, 



from which (see ( |3.35| )) 



Ein(u) = 

v 



(3.41) 



du', (3.42) 



(3.43) 



(3.44) 



(3.45) 



This result, corresponding to the continuous Toda case, has been already obtained 
in [|IIJ]. Consequently, the function T(a, u; | 7 |) defined by ( |3.37| ) can be considered 
as an extended version of the incomplete gamma function T(a,u) = e~ l t a ~ x dt 
nj p. 260). 
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The quantity ( |3.43| ), where T(0,u; | 7 |) can be interpreted as a generalization of 
the exponential-integral function, constitutes an implicit solution of Eq.(l.l). 



4 Linearization and algebraic properties 

Equation (1.1) can be handled within the prolongation scheme ||. In doing so, let 
us consider the prolongation equations for Eq.(l.l): 

yl = F\u,u t -y), y\ — G* (u, u Xj ; y) , (4.1) 
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where i — 1,2, ...N (N arbitrary), and the set of variables y = {y 1 } is the pseudo-po- 
tential @ (j = 1,2, ...M (M arbitrary)). 

The compatibility condition for Eqs. (^4.1|) gives 



7-1 



L i {y)u t + M\u,y), (4.2a) 



\ 7-2 

U 



G* =L\y) [1 + -J u x + P\u; y), (4.2b) 

where M l = M l (u;y), P % = P l (u;y) and L % = L % {y) are defined by 

^±Ml + [P,L\ =eL\ (4.3a) 

(l + -Y [L,M\*=Pi, (4-3b) 



7, 



[M,Pf =0, (4.3c) 



with [P,L]* = P fc g£ - L fe |^, and so on. 

In order to explore the prolongation equations ([4.3| ), for brevity we shall omit the 
index i. 

It is noteworthy the following 



Proposition 3 Let u be a solution of Eg. ( \1.3j ). Then, the function y2 = y2(x,t) 
defined by 

V2x = W (u) sinh y x , y 2t = S(u) cosh y u (4.4) 

7 I u \ 

yix = -u t , yu = 1 + - «i, (4.5) 

7-1 V 7/ 

satisfies the wave equation 

ry \ ( U y- 2 

ym = 1 1 + - V2 XX , (4.6) 

7 \ 7/ 

where W(u) and S(u) obey the linear differential equations of the Bessel type 

7-2 

7 - 1 V ' 7 

and 



W m = ^—U + -\ W, (4.7) 



S»„ = ^ T (l + ^" 2 S + l^(l + ^"s„. (4.8) 
7-1 V 7/ 7 V 7/ 
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Equation (1.6) represents a linearized version of Eq.(1.5). 

To prove this Proposition, let us look for a solution to Eqs. (|4.3|) of the form 



M = W(u)V(y), P = S(u)T(y). 



(4.9) 



Then Eqs. ( }4.3| ) provide 



S = 1 —W u , S u =(l + - 
7 V 7. 



7-2 



(4.10) 



and 



X 2 \ = X 3 , [X x , X 3 ] = X 2 , [X 2 , X 3 ] = 0. 



(4.11) 



withL = Xt, V = X 2 , T = X 3 . 

Equations ( [1.10| ) imply equations ( |4.7|) and ( |4.8|) . By means of the change of 
variable 



z = 2i t 



— ^— ( i + - 

7 - 1 V 7, 



Eq. ( |4.7| ) is transformed into the Bessel equation 



z 2 Z zz + zZ z + \ z 2 - — \ Z = 0, 



V 



where Z = Z(z) is related to W by 



(4.12) 



(4.13) 



W 



u 



1 + 7 * 



On the other hand, Eq. ( |4.8| ) takes the form 



(4.14) 



z Z zz + zZ z 



7-1 

, 7 , 



Z = 0, 



(4.15) 



where 



7-1 



s 



1 + 



u 



7, 



(4.16) 



Z stands for a generic Bessel function of index ±- and ± :L ^-, respectively ( JTB 



pag. 358) 
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The commutation rules (|4.11 ) define the Lie algebra corresponding to the Eu- 
clidean group E 2 in the plane. A realization of ( [4.1 1|) in terms of a two component 
pseudopotential y = (2/1,2/2) is 

X x = d yi , X 2 = sinhj/1^2, X 3 = cosh?/! d y2 . (4.17) 

Therefore, with the help of ( |4.2a| ) and ( |4.2b| ), Eqs. ( |4.1| ) take the form expressed 



by (|4.4| ) and (|4.5|) . Furthermore, by differentiating y 2x with respect to x and y 2t with 
respect to t (see (|4.4|)) and using ( |4.5| ) and ( [4.10| ), one arrives at the wave equation 



Finally, we observe that for 7 = 3, i.e. in the case of the Fermi-Pasta-Ulam 
potential, Eq. ( |4.7|) is led to the Airy equation ([13|], pag. 446). 

W aa = aW y (4.18) 
via the transformation a = (u + 3). 

The prolongation equations ( |4.3| ) offer the possibility of getting a further insight 
into the algebraic structure of Eq.(l.l). 

To this aim, let us write Eqs.( f4.3|) in the form 

2JZlp z + [M,L] = 0, Tz1 z Hm z = eL + [L,P], 

7 7 
[M,P] = 0, (4.19) 

where z = (l + . Then, we look for a solution to Eqs. ( 4. 16 ) such that 

00 00 

M = Y J * k (y)z\ P = J2h(y)z h . (4.20) 

fc=0 fc=0 

In the following, we limit ourselves to characterize mainly the algebraic structure 
of Eq.(l.l) for a particular sequel of values of the parameter 7 (see Proposition 4). In 
this case it turns out that a finite-dimensional Lie algebra is associated with Eq.(l.l). 
This algebra is used to write Eq.(l.l) in a potential form, which allows us to estabilish 
some analogies between the prolongation and the symmetry approaches. In general, 
namely for any value of 7 (provided that 7 ^ 0,1) a systematic analysis of the 
algebraic properties of Eq.(l.l), which is based on the ansatz (5.2), requires further 
efforts. 

Substitution from (|4.20|) into Eq. ( |4.19| ) gives the following constraints between 
the coefficients a k {y) and bk{y). 

k ^^b k + [a fe _!,L] = 0, (4.21a) 

7 

(<2i + 2a 2 z + 3a 3 z 2 + Aa^z 3 + . . .) z~ = 

7 

= eL + [L, b + b lZ + b 2 z 2 + b 3 z 3 + ...], (4.21b) 
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and 



[oq, 60] = 0, 

[a ,6i] + [ai,b ] = 0, 

[00,62] + [01,61] + [a 2 ,6 ] = 0, 

N 

X][ a fc-l;6jV-ft] = 0, 
k=l 

with N arbitrary. 

The following property holds 

Proposition 4 If a 7^ and v = ^ s suc ^ that v . . . 
Eqs. ( \4 ■ 21\ ) imply the finite-dimensional Lie algebra L 

[a ,bo] = [a ,bi\ = [L, 61] = 0, 
[60,61] = e6i, 
[60, L] =eL, 



(4.21c) 



[a ,L] 



u-2 



61. 



-1,0,1,2,3... 

(4.22a) 
(4.22b) 
(4.22c) 

(4.22d) 



for e 7^ 0, and the finite- dimensional Lie algebra £ 

[LA] = [60,6!] = [6 , L) 



h),6 ] 

[a , L] 



h), 61] 
1 

u-2 



61. 



(4.23a) 
(4.23b) 



for e = . 



1, 2, 3, . . .), and 



The proof is straightforward. In fact, under the assumption that v 7^ 
-3, -2, -1, . . . , from ( ^2lTj) we obtain a k = 0, [b k , L] = (fc 
the commutation relation (|4. 22c|) . Furthermore, Eq. (|4.21c|) entails the commutation 
relation ( |4.22d| ) and 62 = 63 = . . . =0. Furthermore, Eq.( f4.22b"D can be deter- 
mined by elaborating Eq. (|4.22c|) via the Jacobi identity applied to [do, [60, L]]. The 
commutation relations ( 4.23|) emerge immediately. 

A matrix representation of £ is 



a 






(e 








60 = 





e 







\0 
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6i 




Then, Eqs. (|4.1|) take the form 



/0 
VO 



(4.24) 





e 






e 



from which yi = Aie e ', y 2 = (Ao — Aix) e e *, and 



7-2 




(4.25a) 



(4.25b) 



V3( 



V3t 



Ai 



7-1 



<9C 



c 



(4.26a) 



(4.26b) 



where ( = Xq — Aix, and Ao, Ai are constants of integration. Here can be 
interpreted as a potential variable. 

Equations ( 4.25a|) and ( 4.25b ) allow us to find, in theory, special explicit solutions 

2 

to Eq.(l.l). For example, let us assume that u t = f3(t)(~<- 2 . Then, Eqs. fl4.26aD and 
(|4.26b|) give rise to the differential equation for (3{t): 



(4.27) 



where b 



(7-2)" 



V. 



27(7-1)' 

Equation ( [4.27] ) can be written as 



V + eV = kV^\ 

1 

) through the position 



(4.28) 



$ + e(3 = V"- 1 . (4.29) 

For e = 0, we may use for the algebra C the same representation (5.6), where 
now bo is the null matrix. In this case Eq. (|4.28| ) coincides with Eq.(3.6), which leads 
to explicit solutions of the original equation (1.5) for special values of the parameter 
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7 (see Sect. 3). We remark that via the change of variable V = e~ w and by a 
suitable choice of the free constant k, Eq. (f4.28|) becomes Eq.(3.4), which arises from 
the symmetry operator X4 in the context of the Lie group theory. This fact suggests 
the existence of a possible link between the prolongation method and the symmetry 
approach. However, our result is indicative only, and many aspects of the problem 
remain to be elucidated. 



5 Comments 

We have dealt with a nonlinear field equation arising as the continuous limit of a 
lattice model containing many cases of physical significance (the harmonic, Toda, 
Fermi- Pasta-Ulam, Coulomb-like lattices, and others). This equation, which is new 
at the best of our knowledge, can be considered as a paradigm for the application 
of different analytical procedures. The addition of a perturbative friction-like term 
has allowed us to check the limits of validity of these methods. These are the sym- 
metry and approximate symmetry approach, and the prolongation technique. The 
joint resort to these methods leads to the discovery of some interesting properties of 
Eq. ( |1 . 1| ) , which have been expounded in the Introduction. 

The spirit of this paper is both of methodological and speculative character. The 
results obtained are a challange to apply the same strategy to investigate a three- 
dimensional version of Eq. ( |1 . 1| ) . This purpose might be important in nonlinear field 
theories based on the deformation of the algebra used in the study of dispersionless 



field equations [19 . 
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